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ABSTRACT 

Ground-based sub-mm and mm-wavelength blank-field surveys have identified more 
than 100 sources, the majority of which are believed to be dusty optically-obscured 
starburst galaxies. Colours derived from various combinations of FIR, submillimetre, 
millimetre, and radio fluxes provide the only currently available means to determine 
the redshift distribution of this new galaxy population. 

In this paper we apply our Monte-Carlo photometric-redshift technique, intro- 
duced in paper I (Hughes et al. 2002), to the multi- wavelength data available for 77 
galaxies selected at 850/im and 1.25 mm. We calculate a probability distribution for 
the redshift of each galaxy, which includes a detailed treatment of the observational 
errors and uncertainties in the evolutionary model. The cumulative redshift distribu- 
tion of the submillimetre galaxy population that we present in this paper, based on 50 
galaxies found in wide-area SCUBA surveys, is asymmetric, and broader than those 
published elsewhere, with a significant high-z tail for some of the evolutionary models 
considered. Approximately 40 to 90 per cent of the sub-mm population is expected 
to have redshifts in the interval 2 < z < 4. Whilst this result is completely consis- 
tent with earlier estimates for the sub-mm galaxy population, we also show that the 
colours of many (<, 50 per cent ) individual sub-mm sources, detected only at 850/im 
with non-detections at other wavelengths, are consistent with those of starburst galax- 
ies that lie at extreme redshifts, z > 4. Spectroscopic confirmation of the redshifts, 
through the detection of rest-frame FIR-mm wavelength molecular transition-lines, 
will ultimately calibrate the accuracy of this technique. We use the redshift probabil- 
ity distribution of HDF850.1 to illustrate the ability of the method to guide the choice 
of possible frequency tunings on the broad-band spectroscopic receivers that equip the 
large aperture single-dish mm and cm-wavelength telescopes. 

Key words: submillimetre, millimetre, cosmology, galaxy evolution, star-formation 



1 INTRODUCTION 

The lack of robust redshift information is one of the out- 
standing problems in understanding the nature of the pre- 
sumed high-redshift galaxies identified in blank-field sub- 
mm and mm wavelength surveys. 

In paper I of this series (Hughes et al. 2002) we de- 
scribed Monte Carlo simulations of 250-500/im surveys from 
the balloon-borne telescope BLAST (Devlin et al. 2001) 
and the SPIRE instrument on the Herschel satellite, as 
well as longer-wavelength ground-based 850/im surveys with 
SCUBA. These simulations included a detailed treatment of 



the observational and evolutionary model-dependent errors. 
The sub-mm colours of galaxies derived from these mock 
surveys demonstrated that it will be possible to derive pho- 
tometric redshifts, with a conservative r.m.s. accuracy of 
Az ~ ±0.5, for thousands of optically-obscured galaxies. 
Thus we can look forward to breaking the redshift deadlock 
which, at the moment, prevents a reliable estimate of the 
evolutionary history of a population that contributes a sig- 
nificant fraction (<; 50 per cent for Ssso^m > 2 mjy) to the 
sub-mm background (Hughes et al. 1998, Smail et al. 2002). 

In this second paper we apply the same Monte Carlo 
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photometric-redshift technique to the existing FIR-radio 
multi- wavelength data for 77 sources first identified in blank- 
field 850pm (SCUBA) and 1.2 mm (MAMBO) surveys. We 
calculate the redshift probability distributions of the indi- 
vidual sub-mm and mm galaxies, taking into account obser- 
vational and model-dependent uncertainties, and thus pro- 
vide a measurement of the cumulative redshift distribution 
for the blank-field sub-mm galaxy population. This redshift 
distribution is completely insensitive to ambiguities in the 
identification of the optical counterparts. 

The structure of the paper is as follows: section 2 gives a 
brief description of the Monte Carlo simulations, the method 
to derive photometric redshifts from the radio-sub-mm-FIR 
colours of galaxies drawn from the mock catalogues, and the 
caveats and advantages of this technique over other redshift 
indicators; section 3 presents individual redshift probabil- 
ity distributions for blank-field sub-mm galaxies. In addi- 
tion, we include a brief discussion of the photometric red- 
shifts for radio and FIR-selected galaxies subsequently de- 
tected at 850/im. Furthermore, we determine the accuracy 
of our method through the comparison of the photometric 
redshifts for those 8 sub-mm sources that have published 
rest-frame optical spectroscopic redshifts. Finally, in sec- 
tion 3 we present the cumulative redshift distribution of 
the blank-field sub-mm population. Section 4 discusses the 
likelihood of a radio detection of previously detected blank- 
field sub-mm sources, based only on the known dispersion 
of the 850pm/1.4 GHz colour in starburst galaxies, radio- 
quiet AGN and ULIRGs, which are believed to be represen- 
tative of the high-z sub-mm galaxy population. Section 5 
describes how the determination of redshift probability dis- 
tributions for individual sub-mm galaxies is motivating the 
efforts to provide spectroscopic confirmation of the redshifts 
at mm and cm wavelengths. Ultimately, these spectroscopic 
redshifts will accurately calibrate the photometric redshifts 
derived from rest-frame FIR-radio-wavelength data. Ap- 
pendix A contains the complete catalogue of spectral en- 
ergy distributions (SEDs) and redshift distributions derived 
in this paper. 



2 PHOTOMETRIC REDSHIFT ESTIMATION 
TECHNIQUE 

The Monte Carlo simulations described in this paper pro- 
duce the redshift probability distribution for an individual 
galaxy. The advantage of this technique, compared to the 
popular maximum-likelihood methods, is that it provides 
more information than a simple estimate of the first and 
second moments of the redshift distributions. A detailed de- 
scription of the technique can be found in paper I; we only 
offer a summary here. Hereafter we will refer to a galaxy de- 
tected in a SCUBA or MAMBO survey as a sub-mm galaxy, 
and a sub-mm galaxy generated in the Monte Carlo simula- 
tions as a mock galaxy. 

We choose an evolutionary model for the 60pm lumi- 
nosity function that fits the observed 850/im number-counts. 
Under this model, and assuming a sub-mm survey of ~10 
deg 2 , we generate a catalogue of 60/im luminosities and red- 
shifts. Randomly selected template SEDs are drawn from a 
library of local starbursts, ULIRGs and AGN, to provide 
FIR-radio fluxes, and hence colours, for this mock cata- 



logue. The fluxes of the mock galaxies include both photo- 
metric and calibration errors, consistent with the quality of 
the observational data for the sub-mm galaxy detected in a 
particular survey. We reject from the catalogue those mock 
galaxies that do not respect the detection thresholds and 
upper-limits of the particular sub-mm galaxy under analy- 
sis. In this paper the redshift probability distribution of a 
sub-mm galaxy is calculated as the normalized distribution 
of the redshifts of the mock galaxies in the reduced cata- 
logue, weighted by the likelihood of identifying the colours 
and fluxes of each mock galaxy with those of the sub-mm 
galaxy in question. Thus, we also take into account as much 
prior information on the sub-mm galaxy population (and 
their uncertainties) as possible. For example, we consider 
uncertainties in the sub-mm galaxy number counts, their 
favoured evolutionary models, the luminosity function of 
dusty galaxies, lensing by a foreground magnifying cluster, 
etc. In the case that the sub-mm galaxy is identified as a 
lensed source, then the mock catalogue is also amplified by 
the appropriate factor, allowing intrinsically fainter objects 
to be introduced into the reduced catalogue, using the same 
observational measurement errors. 

We adopt a flat, A-dominated cosmological model 
(H = 67 kms _1 Mpc _1 , « M = 0.3, Ov = 0.7). To quan- 
tify the sensitivity of the individual redshift distributions 
on the assumed evolutionary history of the sub-mm galaxy 
population, we consider six different models that are able 
to reproduce the observed 850/im number-counts within the 
uncertainties. We will refer to them by the following codes: 

lei: luminosity evolution of the 60pm luminosity function 
<E>[L, 2 = 0] (Saunders et al. 1990) as (1 + z) 3 ' 2 for < z < 

2.2. Constant evolution is imposed between z = 2.2— 10. No 
galaxies beyond z = 10, z > 10] = 0. 

lelL13: same evolutionary description as in lei, but im- 
posing a luminosity cut-off above 1O 13 L . 
Ie2: luminosity evolution of 0] as (1 + z) 3 for < z < 

2.3. Constant evolution is imposed between z — 2.3 — 6. 
<E>[L, z > 6] = 0. 

le2L13: same evolutionary description as in le2, but im- 
posing a luminosity cut-off above 1O 13 L . 

ldel luminosity evolution o/$[L,0] as (1 + z) 1 ' 9 and den- 
sity evolution as (1 + z) 2,5 for < z < 3.5; luminosity and 
density evolution as (1 + z)~ ' for 3.5 < z < 7 

lde2 luminosity evolution o/<3?[L,0] as (1 + z) 3 ' 5 and den- 
sity evolution as (1 + z) 3 for < z < 2.0; no further lu- 
minosity evolution, and density evolution as (1 + z)~ 2 for 
2 < z < 6 

The models used in this analysis provide a wide range of 
the many possible evolutionary scenarios that reproduce the 
sub-mm to mm number counts. These models will allow us 
to test the robustness of the photometric redshifts that we 
will derive in the next section. 

An extended library of SEDs (compared to paper I) 
is used in the Monte Carlo simulations. The library con- 
tains 20 starbursts, ULIRGs and AGN that have been either 
extensively mapped at FIR-radio wavelengths, or galax- 
ies that are sufficiently distant that single-beam photomet- 
ric observations measure the total flux: M82, NGC3227, 
NGC2992, NGC4151, NGC 7469, NGC7771, NGC 1614, 
IZwl, Mkn231, Arp220, Mkn273, NGC 6240, UGC5101, 
IRAS 10214+4724, IRAS 05189-2524, IRAS 08572+3915, 
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IRAS 15250+3609, IRAS 12112+0305, IRAS 14348-1447 
and Cloverleaf. This selection has the important conse- 
quence of reducing any dispersion in the SEDs due to aper- 
ture effects. A continuum SED that includes synchrotron 
and free-free emission in the radio to mm regime, and ther- 
mal emission from dust grains in the mm-to FIR regime, is 
fitted to the observational data (Chapin et al. 2003). The 
luminosities and temperatures of the dust that dominates 
the rest-frame FIR emission are evenly distributed in the 
ranges 9.0 < logL Fm /L < 12.3 and 25 < T/K < 65. 

Half of the galaxies in the SED library used in the 
derivation of photometric redshifts in this paper are in com- 
mon with those considered by Yun & Carilli (2002), here- 
after YC02. The temperature range covered by our sample 
of galaxies (25 — 65K, median 41K) is broader than that of 
YC02 due to the increased fraction of cooler sources, and 
is, for the same reason, colder than those of well-studied 
ULIRGs (46 - 77K, median 60K, Klaas et al. 2001). Despite 
the similar temperatures of our sample to late-type galaxies 
selected at 60/im (25 — 55K, median 35K), drawn from the 
SLUGS catalogue of Dunne et al. (2000), the shapes of the 
SEDs in our library differ significantly from the local galax- 
ies in the SLUGS. This illustrates an important point: it is 
the shapes of the SEDs, and not the derived dust temper- 
atures, that influence the estimated photometric redshifts. 
Given the limited FIR-sub-mm data, there is a degeneracy 
between the dust temperatures, grain emissivity index (f3), 
source-size (£1) and optical depth (t„). This makes it possible 
to fit a single observed SED with a broad range of tempera- 
ture, by tuning the choice of /3, f2 and t„ (e.g. Hughes et al 
1993). 

It remains unclear whether the lower-luminosity, local 
galaxies in the SLUGS are closely related to the blank-field 
SCUBA population which have rest-frame FIR luminosities 
> W 12 Lq, assuming the sub-mm galaxies lie at z > 1. The 
claim that cold SEDs match the number-counts in hydro- 
dynamical cosmological simulations (e.g. Fardal et al. 2002) 
should not be a basis for favouring cold SEDs, since a differ- 
ent evolutionary model with warmer SEDs can fit the counts 
equally well. 

The strongest caveat in any photometric redshift analy- 
sis is the validity of the assumption that the SED templates 
of local starbursts, ULIRGs and AGN are good analogs of 
the high-redshift sub-mm galaxy population. It is not until 
future rest-frame FIR-sub-mm observations from sensitive 
balloon-borne or satellite experiments (paper I) accurately 
measure the variety of SEDs present in the sub-mm galaxy 
population, and spectroscopic redshifts for these are avail- 
able, that it will possible to improve the error estimates 
presented in this paper. 

The technique adopted in this paper differs philosoph- 
ically from those developed in other photometric-redshift 
studies. Our redshift estimates do not rely on the match of 
the observed data to a single template SED, and our errors 
are not drawn from the departure of other local SEDs from 
the adopted standard template, which has become the com- 
mon practice (Carilli & Yun 1999, 2000, hereafter referred to 
as CY99, CY00; Gear et al. 2000; YC02). Instead, we allow 
the test galaxies to be matched against the whole variety 
of template SEDs that can be detected at any given red- 
shift. This point should be emphasized, since although we 
assign our template SEDs at random, not all SED shapes 



(scaled to the same FIR luminosity) are detected with equal 
probability at any given redshift and wavelength. The red- 
shift distributions derived in this paper are often asymmet- 
ric, and broader than those inferred from the single SED 
template analysis. Despite the efforts of some authors to in- 
clude realistic error-bars in their single-template estimates, 
we consider these to have been underestimated: error bars 
derived from 68% of the galaxies in a template catalogue, is 
not equivalent to the 68% dispersion in the fluxes that the 
whole catalogue can create, since some of the outliers con- 
tribute very significantly to the departures from the mean. 

Our redshift estimation method also has the advan- 
tage of including a common-sense evaluation of whether a 
source at a given redshift is simultaneously consistent with 
the observed colours and fluxes. The probability of gener- 
ating a certain luminosity at a certain redshift (i.e. flux) 
is introduced by the luminosity function and the evolution- 
ary model adopted. As shown below, when sub-mm galaxies 
with detections at more than two wavelengths are consid- 
ered, the results are quite insensitive to the differences in 
the evolutionary models adopted to generate the mock cata- 
logue. The consideration of reasonable evolutionary models, 
however, allow plausible colour-based redshift solutions to 
be excluded. 



3 MONTE CARLO BASED PHOTOMETRIC 
REDSHIFTS 

The available data-set for the sub-mm population of galax- 
ies is usually restricted to a low-S/N detection at 850/im 
and an upper-limit at 450/im (from SCUBA surveys), or a 
1.2 mm detection (from MAMBO surveys or an IRAM PdB 
follow-up observation), with, perhaps, radio observations at 
1.4, 5 or 8.6 GHz. Although a number of the SCUBA surveys 
have been conducted in fields previously observed by the In- 
frared Space Observatory (ISO) with the ISOPHOT camera 
at 170 and 90pm, e.g. HDF (Hughes et al. 1998; Borys et 
al. 2002), ELAIS N2 and the Lockman Hole (Scott et al. 
2002; Fox et al. 2002), the ISO measurements are generally 
too insensitive to provide any additional constraint on the 
photometric redshifts. 

The method described in section 2 produces mock- 
catalogues that replicate the sensitivities and calibration ac- 
curacies of the original surveys in which each sub-mm galaxy 
was detected. The colours and fluxes of the mock galaxies 
provide diagnostic diagrams to illustrate the derivation of 
the most probable redshift: colour-colour-redshift (C — C — z) 
diagrams, when more than 3-band detections are available; 
colour-flux-redshift (C — f — z) diagrams, when just two 
bands are available; and fiux-redshift (f — z) diagrams; when 
just one band detection is available. From these diagrams 
we derive the corresponding redshift probability distribu- 
tion (P — z) for each sub-mm galaxy. The distributions for 
a selection of sub-mm galaxies considered in this paper are 
shown in Appendix A. 

3.1 Specific examples: LH850.1 and BCR11 

LH850.1 is the brightest 850pm source identified in the wide- 
area UK 8 mjy SCUBA survey of the Lockman Hole (Scott 
et al. 2002), and has one of the most complete FIR-radio 
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Figure 1. Colour-colour-redshift (C — C — z) plot for LH850.1 
(Lutz et al. 2001, Scott et al. 2002). The flux ratios of the mock 
galaxies, generated under the evolutionary model le2 (section 2), 
are represented as diamonds, and their redshifts are coded in 
colour according to the scale shown to the right of the panel 
(colours can be seen in the electronic version of the paper) . The 
cross represents the measured colours of LH850.1, and the dashed 
box shows the lcr uncertainty in each colour. 



SEDs (Lutz et al. 2001). BCR11, named after source 11 in 
the sub- mm catalogue of Barger, Cowie & Richards (2000), 
is a low S/N, but nonetheless representative, 850/im source 
detected in a follow-up SCUBA photometry survey of the 
micro- Jansky 1.4 GHz radio sources in the Hubble Deep 
Field (Richards 2000). 

Even though BCR11 was not first identified in a blank- 
field sub-mm survey, and hence is not included in the com- 
bined redshift distribution of the sub-mm selected galaxy 
population, it provides a useful source to discuss. The sub- 
mm photometry observations of BCR11 at the position of a 
known radio source (Richards 2000) guarantees that we have 
the correct association of a radio and sub-mm source, and 
hence the correct SED to analyse. BCR11 allows us to confi- 
dently derive the photometric redshift and accuracy that we 
can expect from only the 850/im/1.4 GHz colour, which is 
commonly used to measure redshifts of blank-field sub-mm 
galaxies with follow-up radio observations (CY99, CY00). 

Figs. 1,2 and 3, respectively, show the C — C — z and 
P—z distributions for the multi- wavelength data of LH850.1, 
and the comparison of the observed SED and template SEDs 
(redshifted to the mode of the P — z distribution). Figs. 4,5 



Figure 2. Discrete and cumulative redshift probability distribu- 
tions of LH850.1. The six estimates, plotted with lines of differ- 
ent style, correspond to the six evolutionary models introduced 
in section 2. 
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Figure 3. The observed SED of LH850.1 normalised to the flux 
density at 850/xm is shown as squares and arrows. The arrows indi- 
cate 3cr upper limits. The squares denote detection at a level > 3ct, 
with lcr error bars. The template SEDs (lines) are redshifted to 
z = 2.6, the mode of the redshift probability distribution of most 
models considered (see table 1). The template SEDs at this red- 
shift compatible within 3<r error bars with the SED of LH850.1 
are displayed as darker lines (blue in the electronic version). 
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Figure 4. Colour-flux-redshift (C — f — z) plot for BCR11 (Barger 
et al. 2000). The description of the symbols is as in Fig. 1 



Figure 5. Redshift probability distributions for BCR11. The de- 
scription of the plot is as in Fig. 2 
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and 6 show the same distributions for BCR11. In both ex- 
amples (Figs. 2 and 5) the peak of the redshift probability 
distribution depends little on the details of the evolution- 
ary model used to compute the mock catalogues. In both 
cases, the models that invoke pure luminosity evolution (lei, 
lelL13, le2, le2L13) have more extended wings, since they 
provide a significant population of sources at redshifts be- 
yond 3 that can reproduce the colours of the observed galax- 
ies. On the other hand, the combined luminosity and den- 
sity evolution models (ledl, led2) have a very strong de- 
cline in the density of sources at redshifts beyond 2 and 3.5, 
and thus, the corresponding redshift distributions are much 
more concentrated to values between 2 and 3. Model ledl 
is the most extreme of these two evolutionary models, since 
the very strong negative luminosity and density evolution 
does not allow any source beyond redshift z = 3.5 to be de- 
tected at any reasonable depth. The concentration effect is 
more dramatic on BCR11, where the colour constraints are 
the weakest, and the wings of the redshift distribution are 
strongly dependent on the prior information on the galaxy 
population that has been used. The sharpness of this dis- 
tribution is, thus, due to the prior, and not to the colour 
constraints. In contrast, the sharpness of the distribution of 
LH850.1 is intrinsic to the well determined colours of the 
source, and is almost independent of the prior. 

The comparison of LH850.1 and BCR11 demonstrates 
the improved redshift accuracy that, in general, one can 
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Figure 6. Comparison of the observed SED of BCR11 with those 
of template galaxies, redshifted to z = 2.7 (table 2). The descrip- 
tion of the plot is as in Fig. 3. Note that the selection of SEDs to 
represent BCR11 is different from that in LH850.1. 
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expect for those sources with the greatest combination of 
deep multi- wavelength data. In the case of BCR11, it also 
illustrates the difficulty of constraining individual redshifts 
with observations at 1.4 GHz and 850/im only. The redshift 
probability distributions of both sources have very similar 
modes (2 = 2.6 — 2.8), yet for the less constraining evolution- 
ary models (lei, le2, lelL13, le2L13) there remains 40-50% 
and 25-35% probabilities that BCR11 has a redshift > 3 
and > 4 respectively. The corresponding probabilities for 
LH850.1 are 25-30% and 5% respectively. For models with a 
strong decline in the number of sources with redshift (ledl, 
led2), the high- 2 tails are not produced. The sharpness of 
the derived redshift distributions is due to the lack of bright 
sources at high redshift, and not by the colours of the galaxy. 

Fig. 7 shows the manner in which the redshift distribu- 
tion derived for LH850.1 (using model le2) degrades as dif- 
ferent observational constraints are successively taken into 
account. When the 450, 850/im, 1.2mm and 1.4GHz detec- 
tions are included in the analysis, the addition of upper lim- 
its to the SED makes no difference to the derived redshift 
distribution. It is also apparent that the distribution does 
not change significantly whilst the detection at 450/im is 
used in the analysis, in combination with at least one ad- 
ditional band in the mm-sub-mm and radio regime. How- 
ever, as soon as the 450/im data is excluded, a high redshift 
tail appears. This is due to the large intrinsic scatter in the 
1.4GHz/850/im ratio amongst the template SEDs which, in 
the absence of the 450/im data, allows galaxies over a wide 
range of redshifts to satisfy the colour constraints. The elim- 
ination of the radio data from the analysis allows more low- 
redshift sources to be possible counterparts of the detected 
source, and hence this flattens the resulting redshift distri- 
bution. 

The improvement in constraint on the photometric red- 
shift is extremely important if it is to be used to deter- 
mine the tuning of any broad-band mm-cm spectrograph, 
in an attempt to detect redshifted molecular CO-lines. This 
is discussed further in section 5. Higher frequency radio ob- 
servations at ~ 5 — 15 GHz and deep 450/im observations, 
including sensitive upper-limits, are essential ground-based 
data that can provide the additional diagnostic power to 
constrain redshifts (e.g. LH850.1). 

3.2 Individual redshift distributions 

Before we calculate the cumulative redshift-distribution for 
the 850/tm sources detected in the wide-area (> 200 sq. ar- 
cmin) blank-field SCUBA surveys, which include the Lock- 
man Hole and ELAIS N2 regions (Scott et al. 2002; Fox et 
al. 2002), and the CUDSS fields (Lilly et al. 1999; Eales et 
al. 2000; Webb et al. 2002), it is instructive to discuss some 
of the typical redshift distributions that can be derived from 
data of different quality (i.e. sensitivity and wavelength cov- 
erage) for individual sources. 

Therefore, we compare the redshift distributions of 
wide-area blank-field sub-mm sources with those detected in 
various lensing-cluster surveys (Smail et al. 1997; Smail et 
al. 1999; Bertoldi et al. 2000; Chapman et al. 2002a). We also 
consider radio, FIR and UV-selected objects with follow-up 
sub-mm data (Barger et al. 2000; Scott et al. 2000; Chap- 
man et al. 2000), and two bright sub-mm sources with ex- 
tensive multi-wavelength observations, HDF850.1 (Hughes 



et al. 1998; Downes et al. 1999; Dunlop et al. 2002), and the 
extremely-red object HR10 (Dey et al. 1999). 

Tables 1,2, and 3 summarize the photometric redshifts 
of all the above sources: the modes of their redshift distri- 
butions, 68% and 90% confidence intervals derived from the 
six models studied in this paper, and other relevant infor- 
mation. Appendix A contains the C — C — 2, C — f — z, P — z 
diagrams and SEDs of the galaxies listed in Tables 1 and 2, 
and also a selection of those listed in Table 3 for model le2, 
which is one of the less constraining models considered and 
serves as an illustration of the points made here. 

When we have the most complete multi-wavelength 
data-set (rest-frame radio-sub-mm-FIR, including upper 
limits) then photometric redshifts can be constrained to a 
68% confidence band of width Az ~ ±0.5 (Table 1 and 
Fig.Al). Those sub-mm galaxies for which just two detec- 
tions, and perhaps some shallow upper-limits are available, 
do much worse, with accuracies at a 68% confidence level 
of width Az £ ±1 (Table 2 and Figs. A2). A merit of the 
Monte Carlo photometric-redshift method is that it can de- 
termine redshift probability distributions even in the case 
that the galaxies have just been detected in one band (Ta- 
ble 3 and Fig. A3), the result being very flat redshift distri- 
butions, that are highly dependent on the assumed evolu- 
tionary model. This is especially important because half of 
the sub-mm sources in the blank-field surveys belong to this 
category, although it is worth noting that several of these 
previously reported sub-mm sources have now been found 
to be spurious (Ivison et al. 2002). 

Often, the well constrained distributions are asymmet- 
ric, skewed towards lower redshifts with a high-redshift tail 
(e.g. LH850.1, N2850.2, HR10, ... in Figs. 2 and Al). This is 
a natural consequence of the ability of the data (particularly 
with a deep 450/im detection or upper limit, < 30mJy) to 
strongly reject the possibility of sub-mm galaxies lying at 
2 <, 2, whilst the same observational data cannot reject the 
possibility that the sub-mm galaxies lie at a much higher 
redshift, 2 3> 3. Examples of the power with which the 
450/im detections or upper-limits can reject z <, 2 can be 
seen in the SEDs of LH850.1 and HDF850.1 (Figs. 3 and 
Al), respectively. On the other hand, the bright ISOPHOT 
detections at 170/im (> 120mJy) make a strong case that the 
sources lie at 2 <, 2, for example Nl-40, Nl-64 (Fig. Al) and 
Nl-8 (Fig. A2), whilst the ISOPHOT non-detections arc 
too shallow to offer any additional constraints on the SEDs 
of the majority of the sub-mm population (e.g. LH850.1, 
HDF850.1 in Figs. 3 and Al). 

It is also possible to obtain a bimodal redshift proba- 
bility distribution for some galaxies (e.g. SMMJ00266+1708 
in Fig. Al, BCR33 in Fig. A2), an effect created by two or 
more template SEDs reproducing the same observed colours 
at very different redshifts. 

The figures in appendix A show the insensitivity of the 
850/im/1.4GHz ratio to redshift when the sources are at 
2 £ 2 (a point also made by CY99, CY00). The redshift dis- 
tributions in these cases are usually shallow, with long high- 
redshift tails (e.g. BCR11, BCR13, MMJ154127+6615, in 
Fig. A2), although exceptions also occur (SMM J 16403+4644 
in Fig. 2A). The distributions are even broader in the cases 
when only a sub-mm detection and a radio upper-limit are 
available (Fig. A3). All these redshift distributions are nar- 
rower for those evolutionary models with a strong decline 
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Figure 7. Redshift distributions and compatible template SEDs for LH850.1, when different photometry bands are taken into account. 
All estimates are based on a Monte Carlo produced with the same evolutionary model (le2), The SED models represented in the right-hand 
panels have been redshifted to the mode of distributions represented in the left. Symbol and line codes are as in Fig. 3. 
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beyond z ~ 2 — 3, as can be seen from the measurements in 
Tables 2 and 3, especially for ldel which essentially contains 
no sources at high-z. 

A quick browse through Appendix A illustrates the va- 
riety of template SEDs that reproduce the colours of an 
individual sub-mm galaxy, which lends support to the use 
of a template library in the derivation of the photometric 
redshifts. A few sub-mm sources, however, still cannot be 
reproduced at any redshift with the template SEDs avail- 
able (e.g. LH850.12, SMMJ 16403+46437). A possible expla- 
nation for this problem is the contribution of AGN emission 
to the radio-fluxes at a level not included in the SEDs of 
the template library, or a mis-identification of a sub-mm 
galaxy with a nearby radio-source (see also section 4). A 
wider range of SEDs, including modest fractions of galaxies 
with an excess of radio-emission, could partially solve this 
problem but, on the other hand, it would also introduce an 
extra error term in the redshift computation for the rest 
of the sub-mm galaxies. Other sub-mm galaxies have only 
a small representation of analogs among the mock galax- 
ies (e.g. Nl-40, SMMJ14009+0252), most probably because 
they are rare objects in the Universe, either by their intrinsi- 
cally high luminosity or their low redshift. This latter point 
is illustrated by the rarity of mock galaxies in the corre- 
sponding C — C — zorC — / — z diagrams, assuming any of 
the evolutionary scenarios considered in this paper. 

3.2.1 Accuracy of the photometric-redshift method 

Fig. 8 shows the comparison of photometric and pub- 
lished spectroscopic redshifts for 8 galaxies: in order of in- 
creasing spectroscopic redshift, Nl-40, CUDSS14.18, Nl- 
64, SMMJ02399-0134, HR10, SMMJ14011+0252, 
SMMJ02399-0136, and W-MM11. 

The two galaxies which depart most from the z p hot = 
z Bpcc line are SMMJ02399-0134 and W-MM11. The first 
has a photometric redshift calculated only on the basis of its 
1.4 GHz/850/im colour (Fig. A2), and the latter on a 850/im 
detection with a 450/xm upper limit (Fig. A3). These two 
objects illustrate the large confidence intervals that have to 
be considered when the redshifts are derived from a single 
colour. Specifically, SMMJ02399-0134 illustrates the uncer- 
tainties attached to the large dispersion of 1.4GHz/850//m 
colours implied from a local sample of galaxies. The agree- 
ment between spectroscopic and photometric redshifts for 
the remaining six galaxies (with more than one colour avail- 
able) is remarkably good. The asymmetry of the error bars 
is a reflection of the skewness of their respective P — z distri- 
butions (Figs. Al, A2). These sources illustrate the impor- 
tance of obtaining the description of the whole redshift dis- 
tribution instead of relying on measurements of the median, 
which departs significantly from the most probable redshift 
in skewed distributions. 

A larger sample of spectroscopic redshifts is obviously 
essential to assess the statistical significance of the good- 
ness of this technique, and the accuracy of the error bars 
we derive. For the sources with more than one measured 
colour, these errors are mainly the result of the scatter in 
the shapes of the template SEDs. A larger sample of spec- 
troscopic redshifts would allow us to fine-tune the selection 
of SEDs and priors in order to trim the error bars if, for 
instance, the modes of the distributions do not populate 




spectroscopic z 



Figure 8. Comparison of the photometric-redshift estimates and 
the true redshifts of sub-mm sources with published optical/IR 
spectroscopy. The diamonds represent the modes of the redshift 
distributions for le2. The solid and dashed error bars represent 
68% and 90% confidence intervals, respectively. Sources repre- 
sented in black (in increasing redshift: Nl-40, CUDSS14.18, Nl- 
64, HR10, SMMJ14011+0252, and SMMJ02399-0136) have pho- 
tometric redshifts derived from a combination of colours, and are 
the most precise. Sources represented in grey (in increasing red- 
shift: SMMJ02399-0134, and W-MM11) have photometric red- 
shifts derived from only one colour (or limit), and have shallow 
probability distributions, still compatible with their spectroscopic 
redshifts. Note that the asymmetry of the error-bars is a reflec- 
tion of the asymmetry of the redshift probability distributions, 
which in cases (CUDSS14.18 and HR10) are significantly skewed. 

the envelope of error bars implied from the analysis. The 
six best constrained redshift distributions (shown in black 
in Fig. 7) do show a tendency to cluster more around the 
Zspec = Zphot line than the collective extension of their error 
bars. However, fine-tuning the inputs of the Monte-Carlo at 
this early stage of confirmation of the photometric redshifts 
is premature. 

3.2.2 Notes on individual sources 

LH850.1: Other redshift estimates are 2.95±£:<| based 
on the 1.4GHz/850/nm spectral index technique; and 2.72 ± 
0.37 based on a \ 2 minimization with one template SED 
(YC02) . Our own estimate z = 2.6±o;| (le2) is consistent 
with these measurements. 

LH850.12: No SED can reproduce the observed pho- 
tometry of this source. Ivison et al. (2002) find evidence of 
variability at 1.4GHz, and report associated X ray emission. 
They conclude that LH850.12 might be a radio-loud QSO. 
This kind of objects is actually not represented among our 
SEDs, and thus the redshift we derive is unreliable. Note 
that there are no counterpart mock-galaxies in the vicinity 
of the error box of its C — C — z diagram in Fig. Al. 

HDF850.1: Other redshift estimates are > 2.6, based 
on the 1.4GHz/850/im spectral index technique; and 4.11 ± 
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Table 1. Sources detected in three or more passbands at a > 3<r level. Column 1 gives the name of the source; column 2 the 
references for the SED data; columns 3-8 the mode of the photometric-redshift distributions and 68% confidence intervals for 
different evolutionary models, and in parenthesis below, the corresponding 90% confidence intervals; and column 9 notes for each 
object. The references are coded as follows: (1) Lutz ct al. 2001, (2) Scott et al. 2002, (3) Ivison ct al. 2002, (4) Downes et al. 2001, 
(5) Dunlop et al. 2002; (6) Eales et al. 2000, (7) Smail ct al. 1999, (8) Frayer ct al. 2000, (9) Ivison ct al. 1998, (10) Ivison ct al. 
2000, (11) Dey et al. 1999, (12) Chapman et al. 2002b. (13) Cowie et al. 2002 The A in the notes denotes the amplification factor 
reported in the literature for the lensed sources. 



object 


ref 


-lcl 
^phot 


_lelL13 
phot 


_le2 
phot 


z lc2L13 
phot 


-ldel 
phot 


,lde2 
phot 


notes 


LH850.1 


1,2,3 


o fi+0.4 
Z - D:t 0.6 


2.6±g;t 


2-6±g;t 


2-6±8;| 


2.7±0.3 


9 7+03 








(2.0-4.1) 


(2.0-3.8) 


(2.0-3.7) 


(2.0-3.9) 


(2.0-3.1) 


(2.3-3.5) 




LH850.3 


2,3 


2.7±0.2 


2.7±0.2 


3-0±°;g 


2.5±°;° 


2.9±?6 


3-o±8:§ 








(2.1-3.0) 


(1.8-3.0) 


(2.1-3.0) 


(1.0-3.0) 


(1.0-3.5) 


(2.1-3.0) 




LH850.8 


2,3 


3.8±H 


3.8±J;| 


3.7±J;f 


3.8±J;| 


o c_|_0.0 

j.o=t 6 


30±8;8 








(2.0-6.5) 


(1.6-6.0) 


(2.3-6.0) 


(2.0-5.3) 


(2.6-3.5) 


(2.0-4.3) 




LH850.12 


2,3 


0.6±°; 3 


o-6±8:? 


2+ ' 8 


0-2±°;S 


0.7±0.2 


1.0±0.5 


radio-loud QSO^ 3 - 1 , z p hot is 






(0.1-1.0) 


(0.1-1.0) 


(0.0-1.3) 


(0.0-1.2) 


(0.5-1.0) 


(0.8-1.5) 


unreliable 


N2850.1 


2,3 


2.9±0.6 


2 8+ - 7 


2 - 8± oi 

(2.0-4.1) 


2 8+ - 3 


2 q+0.3 


2 7+0- 3 








(2.0-4.1) 


(2.0-3.7) 


(2.0-5.3) 


(2.4-3.5) 


(2.0-3.5) 




N2850.2 


2,3 


2-3±8;| 


z ' o=c 0.3 


2 14-0.6 


2-3±8:i 


2.3±0.3 


2 6+ 04 








(1.5-3.4) 


(1.5-3.2) 


(1.5-3.5) 


(1.5-3.0) 


(2.0-3.4) 


(2.1-3.0) 




N2850.4 


2,3 


2 - 8± a3 


2.7±0.7 


2 Q+ 1 - 2 


2- 8 ±H 




2 7+0-3 








(1.3-4.5) 


(2.0-4.5) 


(1.2-4.5) 


(1.6-4.5) 


(2.2-3.5) 


(1.5-3.4) 




HDF850.1 


4,5 




4.1±0.6 


4.1±0.6 


4 0+ 07 


q c:+0.0 


4 n+o- 5 


A = 3 






(3.3-5.5) 


(3.1-5.0) 


(3.0-5.1) 


(3.0-5.0) 


(2.7-3.5) 


(3.2-5.0) 




CUDSS14.1 


6 




3.8±g;| 


3.8±g;| 


3.7±0.8 


2.9±0.4 


2 - 8 ±oi 








(2.5-5.2) 


(2.5-5.5) 


(2.0-4.9) 


(2.2-5.0) 


(2.1-3.5) 


(2.0-4.8) 




SMMJ00266+1708 


7,8 


2 7+ 2 - 2 


2 7+ 2 -° 


2 7+1-6 


2.8± 2 ; 2 


2 8+ - 5 


2 q+0.5 


A = 1.6 






(2.1-5.5) 


(2.2-5.5) 


(2.0-5.4) 


(1.0-5.3) 


(2.2-3.5) 


(2.5-5.3) 




SMMJ02399-0136 


8,9 


3.5±° : ° 


3i±8:l 


3.9±0.9 


3-i±8:e 


Q c;_)_0.0 




z S pcc = 2.8, A = 2.4, 






(2.3-5.0) 


(2.0-4.5) 


(2.5-5.0) 


(2.5-4.6) 


(2.2-3.5) 


(2.4-5.0) 


AGN, 3.4cm too low? 






2.8±0.7 


2.9±0.6 


2 - 8± o:I 


2-8±8;5 


2- 8 ±8:l 


2- 8 ±8:! 


disregarding the 3.4cm upper limit 






(2.0-4.5) 


(2.0-4.4) 


(2.0-5.0) 


(2.0-4.3) 


(2.0-3.5) 


(2.0-4.4) 




SMMJ09429+4658 


7,13 




2.8±0.6 


2 - 8=t o:3 


2-s±8:I 


3.0±8;f 


2 7+ ' 3 


A = 1.5 






(2.0-5.6) 


(1.9-4.0) 


(2.0-4.9) 


(1.9-4.0) 


(2.3-3.5) 


(2.0-4.3) 




SMMJ 14009+0252 


10 


4 4+ - 7 


4 1± - 8 


4.1±0.8 


4.1±0.8 




4-i±?:l 


A = 1.5 






(2.7-5.5) 


(3.0-5.1) 


(3.1-5.0) 


(2.5-5.2) 


(2.1-3.5) 


(2.0-5.0) 


AGN contribution? 


SMMJ14011+0252 


8 


2-s±8:i 


2 S+0.5 


2 7+0.6 
Att 07 


2 7+ - 6 


2 7+ - 6 


2 7+ - 3 


Zspec — 2.57, A = 3. 






(2.0-4.5) 


(2.0-4.5) 


(2.0-4.4) 


(2.0-4.5) 


(2.3-3.5) 


(2.0-3.5) 




HR10 


11 


i - 0=t o.o 


i - 0=t o.o 




i-s±8:S 


LO:t 0.0 


i-s±8:S 


Zspcc = 1-44 






(1.5-3.1) 


(1.5-2.5) 


(1.5-2.9) 


(1.5-2.6) 


(1.5-3.1) 


(1.5-2.4) 




Nl-40 


12 


1.2 ± 0.2 


0.7±g;S 


0.6±°;? 


0.6±0.2 


0.7±0.2 


1 0+ ' 5 


Zspcc = 0.45 






(0.6-1.5) 


(0.5-1.5) 


(0.5-1.5) 


(0.3-1.0) 


(0.5-1.0) 


(0.5-1.5) 




Nl-64 


12 


1 2+ 03 


1.2±0.3 


1 2+ - 8 


1 2+ - 3 


i-2±8:i 


1 2+ ' 6 


Zspcc = 0.91 






(0.5-1.8) 


(0.5-1.5) 


(0.5-2.1) 


(0.5-1.8) 


(-) 


(0.7-2.0) 





0.51 based on a x 2 minimization with one template SED 
(YC02). Our own estimate z = 4.1±g;f (le2) is in good 
agreement with these measurements, and is unaffected by 
the inclusion or exclusion of the possible lensing amplifica- 
tion experienced by this source (A w 3, Dunlop et al. 2002) . 

CUDSS14.1: Other published redshift estimates are 
2.01±J;7?, based on the 1.4GHz/850/im spectral index tech- 
nique; and 2.06 ± 0.31 based on a \ 2 minimization with 
one template SED (YC02). Both these estimates are in the 
low-redshift tail of our redshift distribution, which at a 68% 
significance level places this object at z = 3.8±o;g (le2). 

CUDSS14.18: The inclusion of the 450 and 850/im 
fluxes as detections at 2.3 and 2.8<r (instead of upper limits) 
gives a similar distribution to the one shown in Figure A2. 
Our redshift estimate places this object at z — 0.7± i2 with 
a 68% confidence level (le2), in good agreement with its true 
spectroscopic redshift, z sp cc = 0.66. Other redshift estimates 



are 1.12±o;t|, based on the 1.4GHz/850/im spectral index 
technique; and 1.11 ±0.21 based on a x minimization with 
one template SED (YC02). 

SMMJ00266+1708: Other redshift estimates are 
3.49±f;23 based on the 1.4GHz/850/xm spectral index tech- 
nique; and 3.50 ± 0.45 based on a x minimization with one 
template SED (YC02). Our estimate z = 2.7±g;l (le2) is 
consistent with these measurements. The redshift distribu- 
tion is bimodal, due to the range of SEDs used in the redshift 
estimation technique. 

SMMJ02399-0136: The complete SED reported in 
Frayer et al. (2000), Ivison et al. (1998) is inconsistent with 
the SEDs of our template library at any redshift. This is 
reported to be an AGN (Smail et al. 1999). One possible 
explanation for the failure to detect SMMJ02399-0136 at 
3.4cm is a high-frequency break (between 5.5cm and 0.9cm) 
in the synchrotron spectrum. Such sharp spectral features 
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Table 2. Sources detected in two pass-bands at > 3a level. The columns are identical to those described in 
Table 1. The reference codes are as follows: (1) Scott ct al. 2002, (2) Ivison et al. 2002 (3) Barger et al. 2000, 
(4) Eales et al. 2000, (5) Webb ct al. 2002, (6) Smail et al. 1999, (7) Chapman et al. 2002a, (8) Bertoldi et 
al. 2000, (9) Dannerbauer et al. 2002, (10) Chapman et al. 2002b. 
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( c\ n 1 q\ 
(U.U— l.o ) 


fn n n\ 
(U.U— z.U J 


/n n 1 n^ 
(U.U— i.yj 


(U.z— z.U J 




PTTFlt:^ 1 
UlJoDo.l 





9 1 4-0- 4 
2 - 1=t 0.6 


1 84-0-7 


1 74-0- 7 


1 R4-0-5 


1 Q-i_n a 
1 .yzrzu.4 


1 R4-0-4 








(l.U~ 2.0 J 


/n 7 k\ 
t^U. /— z.O ) 


(1.0— z.O ; 


^n n k\ 

(u.y— z.o ; 


(1.0— Z.DJ 


/to 
(l.z— z.OJ 




V_y UUODO.O 


o 


9 14-0-9 


9 04-0 ^ 
z.uzrzu.o 


9 04-0- 8 


9 04-0 ^ 
z.uzrzu.o 


9 K-r-0.4 


1 Q-i_n 4 
1 .yzrzu.4 


Viia4i 1 zlPN-z APN? 






(1 .4— O.U J 


M 4 ^ O 1 * 
l^±.4— O.U J 


/■in q n\ 
(l.U— O.U ) 


(l.Z— O.U ) 


fo n q nA 
(Z.U - O.U j 


/1 9 9 c;A 
( 1 . z — Z . O ) 






o 


1 O4-0- 2 


1 04-0-4 

A - u± o.o 


1 1-L0.8 

1 - i± o.i 


74-°-^ 


1 04-0-4 


1 ^4-0-5 








/n 1 

(U.O— L.O) 


^l.U— Z-U J 


(U.O— Z.U ) 


(U.U— 1.4 J 


(l.U— 1.0 ) 


(U.O— Z.U J 






o 


z.uzrzu.o 


9 fi4_n Q 
z.uzrzu.y 


9 c;4_l n 
z.o zn_L .u 


9 ^4-0 Q 

z.ozrzu.y 


9 7-1-0.8 


9 ^4-0-0 
z.ozrz og 








(1.4— 4.0 ) 


1^1. 0—4.O J 


(1.0—4.0 ) 


(l.U— 4.Z ) 


(i.y— 0.0 ) 


/t n 9 
(l.U— z.OJ 




O UlJooo.10 


o 


9 9-Ufl 7 
Z.2ztU. i 


z.oztU. i 


9 91 0.9 


9 o_i_n 7 
z.ozrzU. / 


z.OztU.O 


9 k_l0.0 








/in q 7\ 
^l.U— 0. / J 


(1.0— 4. U ) 


(1.0— 4. U J 


( ~\ n q 7^ 
(l.U— 0. / J 


(i.y— 0.0 ) 


/TO O c^A 
(l.Z— z.OJ 




oiviivijuzoyy — uio'i 


u 


9 94-0 7 
z.zr^u. 1 


9 ^4-0 R 
z.ozrzu.o 


9 Q4_n r 

z.ozrzu.o 


9 ^4-0 R 
z.ozrzu.o 


9 fi+0.7 


9 mO.O 
2 - 0± 0.8 


4 — 9 c — 1 nc 
SI — Z.O, Z S pec — l.UO 






/t n q 
^i.u— o.O j 


(~\ a a n\ 

(1.4— 4. U ) 


(l.U— 0.0 ) 


(l.U— 0. ( J 


(1. i — o.OJ 


/I 9 9 
(l.z— Z.OJ 




CJMM 71 09 ^fi 4-041 9 

QlViivl J lU^OUTU'ilZ 


7 


9 ^4-0 R 


9 ^4-0 8 
z.ozcu.o 


9 Q4_n r 

z.ozrzu.o 


9 ^4-0 R 
Z . O ZTZ U.O 


9 7-L0.3 
2 "' ± 0.7 


9 c+0.2 
Z - 6± 0.5 


A = 3 2 






yl.O 4.U ) 


M ^ zl 0~\ 

I l.O 4.U^ 


C1 ^—4 0^ 

(l.O 4.U^ 


(l.O 4.U^ 


^Z.U O.OJ 


f"1 7—^ 0"! 




CJA/fTV/T T1 09^74-0,11 


7 


9 04-1-1 


9 04- 1 - 2 


9 14-1-0 


9 04-1-1 
2 " U± 0.5 


z.ozrzu.o 


9 A4-0 A 
Z.4ZTZU.4 


4 — ^7 

+i — O. / 






(1.0—3.5) 


(1.0-3.7) 


(1.0-3.6) 


(1.0-3.6) 


(1.7-3.5) 


(1.6-3.0) 




OlVllVl J 104UO~r4044 


7 


9 94-0 8 
Z.ZzCU.O 


9 94-0 S 
Z.ZzCU.O 


9 14-0-4 


9 94-0 7 
Z.ZziZU. 1 


9 /i _i_0.6 
2 - 4± 0.8 


9 ^4-0 A 
Z.OZiZU.4 


4 — ^ O 
/1 — O.U 






/n ^ q 1 \ 

(^U. O—O.l ) 


(V. 0—0.0 J 


^l.U— o.O j 


(U.O— o.O ) 


(l.Z— o.O J 


/T K Q nA 
(1.0— O.UJ 




QTV/TA/T T1 ft/1 HQ I /Ifi/I/IH 

olvllVl J 104Uo-|-4044U 


7 


Q_Ln 7 
z.oztU. / 


Q_i_n 7 
z.o±U. ( 


9 7_l_n Q 
2. / ±U.o 


Q_i_n 7 
z.oztU. ( 


9 Q-i-0.6 


9 Q4-0-1 

2.9± 7 


/i — O.O 








H 5-4 Tl 


(1.5-4.1) 


(1.5-4.1) 


(2.1-3.5) 


/1 5_3 ii, 

^ _L . O O.l J 




SMM T1 64034-464S7 


7 


2 2+ - 8 


2 3+ - 7 


2 6± ' 4 
z - 0=c l.l 


2 6-t - 4 


2.6±0.5 


2.5zt0.5 


A = 3.2, 1.4GHz high, AGN? 






(1.0-3.3) 


(1.1-3.5) 


(1.0-3.3) 


(1.2-3.5) 


(1.8-3.5) 


(1.5-3.0) 




MMJ154127+6615 


8 






2 8+1- 7 


3.7zt? : 8 


2.7±0.3 


2.8zb0.3 


A = 1 






(2.0-6.0) 


(2.0-5.9) 


(2.0-5.0) 


(2.0-5.0) 


(2.0-3.1) 


(2.0-4.1) 




MMJ154127+6616 


8 






2.8±5;| 


2-Szzbg;! 


3-0±g;t 


9 q+i° 
z.yzrz 04 


A = 1 






(2.0-6.0) 


(2.0-6.0) 


(2.5-5.0) 


(2.1-5.0) 


(2.2-3.5) 


(2.6-5.3) 




MMJ120546-0741.5 


9 


4-7±i:i 




3.8±g;J 


3-8zzb§;§ 


3.0±°,§ 


3.5zbl.O 








(2.5-7.8) 


(2.5-7.0) 


(3.0-6.0) 


(2.5-6.0) 


(2.5-3.4) 


(2.5-6.0) 




MMJ120539-0745.4 


9 


2.8±J;| 


2.8±J;| 


2 - 8=1: o:3 


2-5±J:S 


2.5±0.3 










(2.0-6.3) 


(2.0-6.0) 


(2.0-4.9) 


(2.0-5.3) 


(2.0-3.3) 


(2.0-3.7) 




Nl-008 


10 


OH- - 5 


OH- - 5 


o+o- 4 


0+ 04 


O+O- 4 

u - u=c o.o 


0+ 04 








(0.0-0.5) 


(0.0-0.5) 


(0.0-0.5) 


(0.0-0.5) 


(0.0-0.5) 


(0.0-0.5) 
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Breaking the "Redshift Deadlock" -II 11 



Table 3. Sources detected at 850/im or 1.2mm at a > 3.5c level, and undetected at any other wavelength. 
The columns are identical to those described in Table 1. The reference codes are as follows: (1) Scott et 
al. 2002, (2) Ealcs ct al. 2000, (3) Webb ct al. 2002, (4) Smail et al. 1999, (5) Chapman et al. 2002a, (6) 
Dannerbauer et al. 2002. (7) Ivison ct al. 2002 



object 




ref 


,lcl 
phot 


z 1g1L13 


2 le2 ot 


2 le2L13 


„ldel 
phot 


2 ldc2 


notes 


LH850.2 




1,7 


5.9±i;| 




5.8±5'? 


s-i±?:t 


ci - 0:t 0.5 


A Q+l-2 










(3.5-9.0) 


(3.3-8.5) 


(3.3-6.0) 


(3.3-6.0) 


(2.7-3.5) 


(3.2-6.0) 




LH850.4 




1,7 


5.5±?;° 


5.3± 2 ; 2 




4 -6±J:| 


3 5+°-° 


3-9+-o!4 










(3.0-8.8) 


(3.0-8.5) 


(3.5-6.0) 


(3.3-6.0) 


(2.6-3.5) 


(3.0-6.0) 




LH850.5 




1,7 


5.3±?;§ 


6.1±i;? 


4 7+ 1 - 3 


6.0±? : « 


3 5+0.0 


4.5±1.0 










(3.3-9.0) 


(3.0-8.5) 


(3.4-6.0) 


(3.4-6.0) 


(2.6-3.5) 


(3.0-6.0) 




LH850.13 




1,7 


4.9± 2 t 




5.0±1.0 


5.0±1.0 


3 5+0.0 


4.6±1.2 










(3.0-8.8) 


(3.0-8.4) 


(3.0-6.0) 


(3.3-6.0) 


(2.6-3.5) 


(2.7-6.0) 




N2850.3 




1,7 




A fi+ 2 -9 
4 ' D=C 0.6 


5.8±°/ 2 


4.8±0.9 


3-5±og 


4.3±1.0 










(3.5-9.0) 


(3.5-8.8) 


(3.6-6.0) 


(3.6-6.0) 


(2.7-3.5) 


(3.2-6.0) 




N2850.5 




1,7 


5.0±2.0 


4.7±?i 




4-5±J;| 


3.5+. ,; , 


2 -5± 2 ;g 










(2.5-8.6) 


(2.7-8.5) 


(2.9-6.0) 


(2.9-6.0) 


(2.5-3.5) 


(3.0-6.0) 




N2850.6 




1,7 


5.0±2.0 


4.8±l;§ 


4.8±?;| 


4.8±1.2 




2 -5±o:l 










(2.5-8.7) 


(2.5-8.3) 


(3.0-6.0) 


(2.5-6.0) 


(2.5-3.5) 


(2.0-5.6) 




N2850.8 




1,7 




4.7± 2 -° 


4.5±};° 


4.8±1.2 


3.0±g;5 


2.5±?;g 










(2.1-8.5) 


(2.0-8.4) 


(2.6-6.0) 


(2.6-6.0) 


(2.3-3.5) 


(2.5-6.0) 




N2850.9 




1,7 


5.2±f° 


5-0±?;g 


5.9±°;| 


5.0±1.0 


3 5+°-° 


3 7+1- 5 










(3.0-9.0) 


(2.8-8.5) 


(3.2-6.0) 


(3.1-6.0) 


(2.6-3.5) 


(2.5-6.0) 




N2850.10 




1,7 


4.7± 2 ;| 


4.5± 2 ;5 


4.7±?;8 


4.8±1.2 




2.5±°J 










(2.1-8.5) 


(2.1-8.5) 


(2.6-6.0) 


(2.6-6.0) 


(2.3-3.5) 


(1.5-4.8) 




N2850.ll 




1,7 


5.5±1.9 




5.0±1.0 


4 8 1,2 


3 - 5± o!? 


4.1±1.2 










(3.0-8.9) 


(3.0-8.7) 


(3.2-6.0) 


(3.2-6.0) 


(2.6-3.5) 


(2.6-6.0) 




CUDSS14.2 




2 


5.2± 2 ; 4 


5.0± 2 ;5 


4 7± - 8 
' ^1.0 


5 - 0± l'5 


3.5±0;° 


2 - 5± i:, 2 , 










(2.5-8.7) 


(2.5-8.8) 


(2.9-6.0) 


(2.8-6.0) 


(2.5-3.5) 


(1.5-5.1) 




CUDSS14.4 




2 


5.0±f| 


5-i±?:§ 


4 6+ 1 - 4 


4 7± 1 ' 3 


3 5± 












(2.5-8.8) 


(2.5-8.8) 


(2.9-6.0) 


(2.8-6.0) 


(2.4-3.5) 


(1.5-5.2) 




CUDSS14.5 




2 


4.8±?;| 


4.9± 2 '1 


4 8+ 1 - 2 


5.3±?:5 


3 5+0-0 


n 5+1.9 










(2.8-9.0) 


(2.8-9.0) 


(3.0-6.0) 


(3.0-6.0) 


(2.5-3.5) 


(2.0-5.7) 




CUDSS14.6 




2 


4.7±?-| 


4.7±?;° 


5.2±?;| 


5.2±5'| 


q 5+0-0 












(2.5-8.8) 


(2.5-8.8) 


(2.9-6.0) 


(2.9-6.0) 


(2.5-3.5) 


(1.5-5.3) 




CUDSS14.8 




2 


4.7±f 7 


4.7± 2 ;° 


4 5+ 1 - 3 


4 7+ 1,9 


3.5± ; g 


2.5±?;8 










(2.5-9.0) 


(2.5-9.0) 


(2.6-6.0) 


(2.6-6.0) 


(2.3-3.5) 


(1.5-5.0) 




CUDSS14.10 




2 


4.7± 2 -^ 


4.8± 2 ;| 


3 9+ 2 -0 


4.3±1.2 


2.9±0.5 


2.5±?;5 










(2.5-9.0) 


(2.5-9.0) 


(2.6-6.0) 


(2.6-6.0) 


(2.3-3.5) 


(1.5-4.8) 




CUDSS3.2 




3 


4.5±i;l 




3.5±};J 


3.9±1.4 


2.9±0.5 


2-5±0;0 










(2.0-8.5) 


(1.9-8.5) 


(2.0-5.8) 


(2.2-6.0) 


(2.2-3.5) 


(1.0-4.0) 




CUDSS3.3 




3 


3-8±U 


3.7± 2 4 


4-0±};§ 


4.0±i;| 


2.9±0.4 












(2.0-8.4) 


(2.0-8.4) 


(2.2-6.0) 


(2.0-5.8) 


(2.1-3.5) 


(1.0-4.0) 




CUDSS3.4 




3 


3.2±?;1 


3.7± 2 -^ 


4.3±?;1 


3.9±1.4 


2.9±0.4 


2.5±?;|S 










(2.0-8.5) 


(2.0-8.3) 


(2.3-6.0) 


(2.2-6.0) 


(2.2-3.5) 


(1.0-4.0) 




CUDSS3.5 




3 


3.7±f 7 


3.7± 2 -| 


3.8±};5 


3.9±1.4 


2.9±0.5 


2.5±?;|] 










(1.8-8.5) 


(1.8-8.5) 


(2.0-5.9) 


(2.0-5.8) 


(2.1-3.5) 


(1.0-3.9) 




CUDSS3.6 




3 


3.7±ff 


4 0+ 2 - 4 


3.7±;-5 


3.9±1.4 


2.9±0.4 


2.5±?;0 










(2.0-8.7) 


(2.0-8.8) 


(2.0-5.8) 


(2.0-5.8) 


(2.2-3.5) 


(1.0-3.9) 




CUDSS3.7 




3 


3.7±?;« 


3.7±?;« 


4.0±l;| 


3.7±1.3 


2-s±8:i 


2.5±? : |] 










(2.0-8.4) 


(2.0-8.1) 


(2.3-6.0) 


(2.0-5.7) 


(2.0-3.5) 


(1.4-4.5) 




CUDSS3.11 




3 


3.8±?;0 


4-5± 2 :£ 


3.8±i; 7 


3 7+1-8 


2.9±0.4 


2.5±?;° 










(1.9—8.5) 


(2 0—8 5) 


(2 3—6 0) 


(2.0—5.7) 


(2.2—3.5) 


(1 0—4 1) 




CUDSS3.12 




3 


3.7± 2 -6 


3.2± 2 ;l 


3.5±};g 


3.9±1.3 


2.9±0.5 


2.5±°;° 










(1.9-8.5) 


(1.9-8.5) 


(2.0-5.8) 


(2.0-5.8) 


(2.2-3.5) 


(1.0-3.9) 




CUDSS3.13 




3 


4.0± 2 ;3 


3-2±?;l 


3.5±};g 


3.9±1.4 


2.9±0.5 


2.5±?;[J 










(1.8-8.5) 


(1.7-8.5) 


(2.0-5.8) 


(2.1-6.0) 


(2.1-3.5) 


(1.0-4.0) 




CUDSS3.16 




3 


3-s±?:§ 


3.2± 2 ; 2 


3.7±i;5 


3.5±J;9 


2.9±0.4 


2.5±°;° 










(1.6-8.5) 


(1.6-8.5) 


(2.1-6.0) 


(2.0-5.9) 


(2.2-3.5) 


(1.0-3.8) 




SMMJ02400- 


-0134 


4 


5.0± 2 -| 


5.2± 2 ;3 




5.2±?;8 


3 5+0-0 


4.0±°;g 


A > 1.9 (2.5?) 








(3.0-9.1) 


(3.0-9.1) 


(3.1-6.0) 


(3.1-6.0) 


(2.5-3.5) 


(2.0-5.5) 




SMMJ04431+0210 (N4) 


4 


4.5±?;| 


4-6± 2 :| 


4 - oi i.o 


A 9+1-6 


3 5+0-0 


2.4±1.0 


A > 1. 








(2.5-8.8) 


(2.5-8.7) 


(2.7-6.0) 


(2.7-6.0) 


(2.5-3.5) 


(1.5-4.9) 




SMMJ04541- 


-0302 


5 


4.7± 2 -0 


4.0± 2 | 


4.7±1.2 


4.7±?;| 


3.0±g;§ 


2.5±?;g 


A = 2.6 








(2.2-8.5) 


(2.3-8.5) 


(2.6-6.0) 


(2.6-6.0) 


(2.4-3.5) 


(1.5-4.7) 
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Table 3. (cont.) 



object 


ref 


lei 
phot 


lelL13 
phot 


z lc2 
phot 


le2L13 
phot 


ldel 
phot 


lde2 
phot 


notes 


SMMJ04542— 0301 


5 


4.7±2.2 




* n > 1 7 


4 - 5 ±r.o 


Q Q_1_0 2 


i-i - i n Q 

2.5±°;9 


A = 4.5 






(2.5-8.9) 


(2.1-8.5) 


(2.7-6.0) 


(2.6-6.0) 


(2.4-3.5) 


(1.5-4.9) 




SMMJ04543+0256 


5 


3 - 9± ?:I 


3.8±f;| 


j ^ i 1 r -\ 

4.0±j J 


4.0±J'5 > 


2.9±0.4 


2 - 5± ?:o 


A = 3.1 






(2.0-8.4) 


(2.0-8.5) 


(2.4-6.0) 


(2.0-6.0) 


(2.0-3.5) 


(1.5-4.5) 




SMMJ04543+0257 


5 


A i-r i 9 9 


4 - 5± ?:s 


4.2±J"J 


a i- i i n 
4.5±J;0 


3.5±j];{j 


2 - 5 ±r.o 


A = 3.4 






(2.0-8.2) 


(2.5-8.3) 


(_2.d-d.UJ 


(2.D-D.0) 


(2.3—3.5) 


(1.5—4.8) 




SMMJ10237+0412 


5 


4.0± 2 ;§ 


3.9±f;0 


4.5±°;£ 


j w i n s 
4 7+ 


2.9±0.4 


2.5±°;° 


A = 3.5 






(2.0-8.7) 


(2.0-8.7) 


(2.0-6.0) 


(2.6-6.0) 


(2.2-3.5) 


(1.0-4.1) 




SMMJ14573+2220 


5 


4.0±?-I 

1.5 


4.2±?-§ 
1.7 


4.0±H 

1.0 


4.2±H 
-^0.8 


2.9±0.4 


2.5±?i 
1.0 


A = 2.8 






(2.0-8.5) 


(2.0-8.6) 


(2.0-6.0) 


(2.5-6.0) 


(2.2-3.5) 


(1.5-4.7) 




SMMJ17223+3207 


5 


4.0±f -| 


4.2±ff 


3.9±};g 


4 -o±ii 


3-0±8;l 


2.5±0-i 


A = 2.8 






(2.0-8.4) 


(2.0-8.4) 


(2.4-6.0) 


(2.5-6.0) 


(2.2-3.5) 


(1.5-4.6) 




SMMJ22471-0206 


4 


4 8+ 2 - 4 


4 8+ 2 - 3 


4 7+ 1 - 3 


4.9±1.1 


3 4+0-1 


2 9+ 1 - 3 
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are not taken into account in the fitting of the radio spec- 
trum. We note however that excluding the 3.4cm upper-limit 
from the modelling of the multi-wavelength data provides a 
redshift estimate z = 2.8± ;t (le2) that is in good agreement 
with the true redshift of SMMJ02399-0136: z = 2.80. Other 
redshift estimates are 1.65±g;l§ based on the 1.4GHz/850/tm 
spectral index technique; and 2.83 ±0.38 based on a % 2 min- 
imization with one template SED (YC02). 

SMMJ09429+4658 (H5): Other redshift estimates 
are > 3.6 based on the 1.4GHz/850/im spectral index tech- 
nique; and 3.86 ± 0.49 based on a \ 2 minimization with one 
template SED (YC02). Our estimates actually place this 
source at a slightly lower redshift z = 2.4±q;| (le2). Given 
the sub-mm detection, the low flux density at 1.4 GHz ex- 
cludes almost all of the template SEDs in our library. 

SMMJ14009+0252: Other redshift estimates are 
1.27±o;54 based on the 1.4GHz/850/im spectral index tech- 
nique; and 1.30±0.23 based a on a \ 2 minimization with one 
template SED (YC02). However large discrepancies (~ 3a) 
exist between the template SED of YC02 and the observed 
450/im and 1.4GHz fluxes, suggesting the fit will have a high 
reduced-x 2 , and that their errors (inferred for good matches) 
might be under-estimated. The flattening of the spectrum 
at 450/xm, relative to the Raleigh- Jeans tail defined by the 
850/im and 1.35mm fluxes, supports a higher redshift esti- 
mate. Our own fit indicates that this object has a redshift 
z ~ 4.1 ±0.8 (le2). The discrepancies between different red- 
shift estimators might be attributed to the possible AGN 
nature of this object (Ivison et al. 2000) . 

SMMJ14011+0252: Other redshift estimates are 
2.53±o;fi2 based on the 1.4GHz/850/im spectral index tech- 
nique; and 2.73 ± 0.37 based on a x 2 minimization with one 
template SED (YC02). Our own estimate is z = 2.7±o;? 
(Ie2) in good agreement with these estimates, and with the 
spectroscopic redshift of the optical counterpart z spoc = 
2.57. 

HR10: Other redshift estimates are 1.49±8;Ii based on 
the 1.4GHz/850/im spectral index technique; and 1.75±0.28 
based on a \ 2 minimization with one template SED (YC02). 



Our own estimate z — 1.5±o;o Qe2) is in good agreement 
with their measurements, and with the true spectroscopic 
redshift of this source z spoc — 1.44. 

SMMJ02399-0134: Other redshift estimates are 
1.17±o;45 based on the 1.4GHz/850/im spectral index tech- 
nique; and 0.93 ± 0.19 based on a \ 2 minimization with one 
template SED (YC02). The spectroscopic redshift of this 
source (z S pcc = 1.06) is inconsistent with our best estimate 
of the redshift, which within a 68% confidence interval places 
it at z = 2.3 ±0.8 (le2). Redshifts in the range 1.0 < z < 1.5 
have a 10% integrated probability. 

SMMJ16403±46437: The detection at 1.4GHz 
(593/iJy) is high when compared with the 850/rni level, and 
there is no acceptable correspondence with any of the SEDs 
in our template library. This, plausibly, might be a misiden- 
tification or a radio-quiet AGN, with enhanced radio lumi- 
nosity compared to the rest-frame FIR luminosity. 

SMMJ16403±46440: Treating the 1.35mm observa- 
tion as a detection (2.2<r) provides no additional constraint 
and gives a similar redshift distribution to the one derived 
from the 450/im/850/im colour and the 1.4GHz and 1.35mm 
upper limits, shown in Figure A2. 



3.3 Cumulative redshift distribution of the 
sub-mm galaxy population 

It is straightforward to calculate the cumulative redshift dis- 
tribution for the sub-mm galaxy population as the coaddi- 
tion of the individual probability distributions. It should be 
noted that the calculation carried out in this paper has the 
advantage of including the whole redshift probability dis- 
tribution of each individual galaxy, and not just the mean 
values of the distributions. We have included in this cal- 
culation of the cumulative distribution only those sub-mm 
sources identified in wide-area blank-field surveys (UK 8mJy 
survey and the CUDSS) at a > 3.5cr level: 50 sources. We 
do not include sub-mm galaxies identified in surveys car- 
ried out towards lensing clusters, or targeted observations 
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Figure 9. Combined redshift distribution for the population of 
sub-mm galaxies in the 8mjy SCUBA survey with one or more 
measured colours (the 14 galaxies contained in Tables 1 & 2). The 
different lines represent different evolutionary models assumed for 
the sub-mm galaxy population (section 2). Upper panel: discrete 
bin probabilities. Lower panel: cumulative distribution. 



of other populations of galaxies (radio, ISO, Lyman-break 
galaxies or extremely red objects), or catalogs that have not 
been fully published (MAMBO surveys), since they can dis- 
tort and bias the characteristic redshift distribution of the 
sub-mm galaxy population. 

Half of sources in complete flux-limited sub-mm sam- 
ples, such as the UK 8mJy and CUDSS surveys, have a 
single 850pim detection with one, or more, additional upper- 
limits at other wavelengths: 28/50 sources detected (> 3.5<r) 
at 850/im belong to this category. Since the redshift distribu- 
tions of these sources are the most dependent on the priors 
used, the interpretation of the combined redshift distribu- 
tions should be based on the range of results given by the 
different evolutionary models. 

Fig. 9 shows the cumulative distribution and the com- 
bined redshift distribution of the 14 galaxies in the UK 8mJy 
survey with redshift distributions derived from SEDs which 
contain at least 2-band detections (objects in Tables 1 and 
2). This figure indicates that 65-90 per cent of the sub- 
mm galaxies could have redshifts in the range 2 <5 z <J 4. 
The results are consistent with the expectation that sub- 
mm galaxies detected in the SCUBA (and also MAMBO) 
surveys represent a high-z population (Dunlop 2001, Blain 
et al. 2002). 



Figure 10. Combined redshift distribution for the population of 
all 50 sub-mm galaxies detected in wide-area surveys at 850/^m 
(> 3.5ct) and listed in Tables 1,2, and 3. The different lines repre- 
sent different evolutionary models assumed for the sub-mm galaxy 
population (see text). Beware that the extended tail towards high- 
redshifts is due to sources which are just detected at 850/tm. The 
upper limits at other wavelengths do not usually help to constrain 
their redshifts between 2 < z < 10, and hence their probability 
distributions are very flat in those regimes. Upper panel: discrete 
bin probabilities. Lower panel: cumulative distribution. 



The sources in the CUDSS sample with one or more 
measured colours (9 sources) appear to have a significantly 
different lower redshift distribution, with 40 per cent of the 
galaxies at z <J 2 and approximately 50 per cent of the 
galaxies at 2 <, z <, 4. The high redshift tail of the CUDSS 
sample is similar to that shown by the 8m Jy survey with ~10 
per cent probability at z > 4. The two galaxies responsible 
for the low-redshift difference are CUDSS3.1 and 3.9, both 
with relatively high 450/im flux densities S^o^m > 77 mjy, 
compared with the 30-40 mjy detection levels for the UK 
8mJy survey sources. Taking the data for the CUDSS sample 
at face- value suggests that it includes an over-density of low- 
z galaxies at z Ss 2, compared to the galaxies in the UK 8m Jy 
survey, and vice-versa. This discrepancy, however, could also 
be explained by a difference in the flux-calibration of the 
450/im data in the two surveys. 

Fig. 10 shows the cumulative redshift distributions of 
the sub-mm galaxy population for the 50 galaxies in the 
8mJy and CUDSS surveys detected at a > 3.5<r level. This 
extended sample includes the 27 sub-mm galaxies with a 
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Figure 11. Distribution of intrinsic fluxes at 1.4GHz correspond- 
ing to a population of galaxies discovered at 850/tm, at a level 
> 8m Jy. Depending on the evolutionary model adopted, 80% of 
the sources would have a flux in excess of 2 to 5/iJy, and hence 
in order to conduct a parallel 1.4GHz imaging follow-up program 
with 80% completeness, the la depth of the images should be 1 
to 2/iJy. 



single detection at 850/im and one, or more, upper limits 
at other wavelengths (Table 3). The combined redshift dis- 
tributions, for the different evolutionary models under con- 
sideration, imply that the majority of the sub-mm galax- 
ies, approximately 50 - 90 per cent, lie between z — 2 — 4. 
The remainder of the galaxies, ~ 10 per cent, lie at z < 2. 
Fig. 10 also shows that 50 per cent of the galaxies have 
colours that are consistent with z > 4. The distribution of 
the high-redshift tail is very sensitive to the choice of evo- 
lutionary model, since the individual redshift distributions 
for the 27 sources with only a single detection are neces- 
sarily very flat over the redshift range z J> 2 (see Fig. A3). 
For these galaxies, the constraint on their redshift distribu- 
tions comes mainly from the flux detected at 850/im and the 
ability to reproduce it with the evolution of the 60/mi local 
luminosity function. More extensive photometry, especially 
including shorter sub-mm wavelengths, is crucial to better 
constrain these flat distributions, as illustrated by Fig. 7 (see 
also Fig. 6 in paper I). 



4 THE LIKELIHOOD OF MIS-IDENTIFYING 
THE RADIO COUNTERPARTS OF HIGH-Z 
SUB-MM GALAXIES 

The strength of the correlation between the FIR and radio 
luminosities in starburst galaxies (Helou et al. 1985, Con- 
don 1992) has supported the expectation that all high-z 
star forming sub-mm galaxies should also have a detectable 
radio-counterpart, provided the radio observations are deep 
enough. 

Although originally encouraged by this assumption, 
some of the deepest VLA searches for radio emission as- 
sociated with SCUBA and MAMBO sources have met with 
varying degrees of success. 

The deep VLA surveys of the Lockman Hole and ELAIS 



N2 fields (Ivison et al. 2002) have a la noise limit > 5 — 9/iJy 
at 1.4 GHz. These surveys have identified 11/16 and 8/14 
sub-mm sources at 1.4GHz, respectively, with peak radio 
fluxes above 4<r and integrated fluxes above 3a, after reject- 
ing 5 and 1 sub-mm sources from the parent samples, re- 
spectively, because of their high sub-mm noise levels. From 
these, there are 9/16 and 5/14 sources in the Lockman Hole 
and ELAIS N2 surveys, respectively, with a final S/N>3 
when one considers the additional errors due to calibration 
and fitting a gaussian to the peak of the radio emission. 
If one limits the analysis to sources > 8 mjy at 850/im, 
then 4/10 and 4/8 sources have been simultaneously de- 
tected with S/N> 3 at 850/im and 1.4GHz. This is in clear 
contrast to the claims of previous authors (Fox et al. 2002) 
that bright sub-mm sources (> 8mJy at 850/im) should be 
detected at flux densities of ~ 100/iJy at 1.4 GHz. 

In order to understand the properties of the radio asso- 
ciations, we have calculated the probability of detection at 
1.4GHz for a mock 8m Jy survey, following the same method- 
ology as in paper I, but using the extended library of SEDs 
(section 2) . The intrinsic 850/tm fluxes of the galaxies are 
convolved with la ~ 2.6mJy measurement errors, and 9 per 
cent calibration uncertainty. Fig. 11 shows the distribution of 
intrinsic 1.4GHz fluxes for the detected sources (including 
errors) above the 8mJy threshold at 850/im. 

We conclude that future surveys at 1.4 GHz must reach 
a 3a sensitivity of 2-5/iJy if radio observations are to pro- 
vide counterparts to 80% of the high-z sub-mm galaxies with 
850/im flux densities > 8 mjy, according to the redshift dis- 
tributions assumed in the alternative evolutionary models. 
Thus, current radio observations have insufficient sensitivity 
to detect the bulk of the sub-mm population with adequate 
S/N, and we must wait for the next generation of new and 
upgraded facilities (eg. e-VLA, e-MERLIN, SKA). 

In the case of the UK 8mJy survey, Fig. 10 shows that 
40-50% (Lockman Hole) and 50-60% (ELAIS N2) of the 
sources should have an intrinsic 1.4GHz flux lower than the 
3a threshold, in the absence of any radio-loud AGN emis- 
sion. These levels agree well with the number of radio coun- 
terparts to > 8mJy sources (Scott et al. 2002, Ivison et al. 
2002), bearing in mind the actual variation in sensitivity 
across the SCUBA wide-area survey fields. 

This potential for misidentifying the correct radio 
sources associated with blank-field sub-mm galaxies will in- 
evitably lead to incorrect estimates of photometric redshifts 
if they are derived in part from the 1.4 GHz fluxes. In the cal- 
culation of the cumulative redshift distribution (section 3.3) 
we have treated all of the published associations of 1.4 GHz 
and 850/im sources as the correct radio counterparts to the 
sub-mm galaxies. 

Despite the rapid acceptance of the radio-sub-mm flux 
density ratio, 0:345 (CY99, CY00), as the preferred method 
with which to estimate the redshift distribution of the sub- 
mm galaxy population, this diagnostic ratio loses its power 
to discriminate between photometric redshifts at z > 2 
(CY99, CY00). We have already demonstrated in paper I 
(Fig. 11) the detrimental effect of the large dispersion in 
radio-sub-mm colours of the local template galaxies on the 
accuracy of the derived photometric redshifts of the sub-mm 
blank field sources. Furthermore, the existing radio interfer- 
ometers will remain limited in their ability to map the large 
areas ~ 1-40 sq. deg. that will be surveyed by the next 
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Figure 12. Redshift probability distributions for HDF850.1. The 
different line types and colours correspond to the different evo- 
lutionary models considered. Model ledl is not plotted, since it 
excludes sources at z > 3.5 due to the strong prior decline of 
the density and luminosity of sources. The unshaded area corre- 
sponds to the redshift range where the 12 CO (J = 1-0) line lies 
in the observable K-band window (18.5 - 26.0 GHz), accessible 
with the 100m GBT. There is an 86-90% probability of detecting 
the J=l-0 CO-line in this GBT window 

generation of sub-mm and FIR experiments (e.g. BLAST, 
Herschel, SIRTF) down to the necessary sensitivity levels 
(3<r ~ 2 — 5/iJy at 1.4GHz), making the use of a^ 4 | an im- 
practical measure of redshift for the majority of the popula- 
tion. Deep interferometric radio observations, however, re- 
main a practical way to pinpoint the source of the sub-mm 
emission with sufficient positional accuracy to investigate 
galaxy morphologies and exploit optical/IR spectroscopic 
line diagnostics. 

This limitation of the radio-sub-mm technique 
prompted the design and construction of a Balloon-borne 
Large Aperture Submillimetre Telescope, BLAST (Devlin et 
al. 2001). The primary scientific goal of BLAST is to break 
the current "redshift deadlock" by providing photometric 
redshifts with sufficient accuracy, without the prior neces- 
sity for secure optical, IR or radio counterparts. BLAST will 
use rest-frame sub-mm and FIR colours to derive redshifts 
with an accuracy Az ~ ±0.5 for > 1000 galaxies in the red- 
shift range < z < 6 detected in wide-area (> 0.5 — lOdeg 2 ) 
surveys at 250, 350 and 500^m (Hughes et al. 2002). 

This redshift accuracy is sufficient to guide the choice of 
possible front-end spectrographs, and their receiver tunings, 
on large single-dish ground-based mm and cm telescopes 
(e.g. 100-m GBT-Green Bank Telescope and 50-m LMT- 
Large Millimetre Telescope) in the effort to detect redshifted 
CO emission-lines. We conclude this paper with the descrip- 
tion of one application of this technique. 



5 BROAD-BAND MM-CM SPECTROSCOPY 

Until recently the brightest sub-mm source in the Hubble 
Deep Field, HDF850.1 (Hughes et al. 1998), has avoided 
the detection of its X-ray, optical, IR or radio counterpart, 
despite an accurate interferometric position at 236 GHz 



(Downes et al. 1999) and the unprecedented depth of 
the multi-wavelength HDF surveys (Williams et al. 1996; 
Thompson et al. 1999; Richards 1999; Hornschmeier et al. 
2000). It is only after the accurate subtraction of the ellip- 
tical galaxy, 3—586, at z ~ 1, however, that a faint lensed 
IR counterpart (K = 23.5, H - K = 1.5, I - K > 5.1) to 
HDF850.1 has been discovered (Dunlop et al. 2002). 

This IR detection now gives greater confidence to the 
likelihood that the coincident radio emission, with a flux 
density of 16±4/iJy at 1.4 GHz, from the combined data of 
MERLIN and the VLA, is genuinely associated with the 
sub-mmm source (Dunlop et al. 2002). Furthermore this 
strengthens the suggestion that a marginal 3.3<r detection 
at 8.4 GHz (VLA 3651+1226, Richards, priv. comm.) that 
may also associated with HDF850.1. 

The redshift probability distribution for HDF850.1, 
based on the multi-wavelength 170/im-1.4 GHz data, indi- 
cates z is comprised between 3.3 and 5.4 with a 86% to 90% 
confidence depending on the adopted prior (Fig. 12). This 
high photometric redshift supports the earlier decision to 
reject the z ~ 1 elliptical galaxy 3—586, as the optical coun- 
terpart. Unfortunately the IR counterpart, HDF850.1K, is 
too faint to attempt IR spectroscopy with even the 10-m 
class of telescopes. 

The most efficient means to determine the spectroscopic 
redshifts of optically-obscured sub-mm galaxies is the de- 
tection of millimetre and centimetre wavelength molecular 
lines using large single-dish telescopes (e.g. GBT or LMT) 
with broad-band receivers (Townsend et al. 2001, Carilli & 
Blain 2002) . In this way, it is possible to eliminate the neces- 
sity for an optical/IR counterpart, or a position accurate to 
< 3arcsec. Given the predicted accuracy of the photomet- 
ric redshift for HDF850.1, and taking all of the uncertainties 
discussed in section 2 into account, we demonstrate in Fig. 12 
that the J=l-0 CO line in HDF850.1 has a 86-90% probabil- 
ity of detection in the K-band receiver (18.0-26.5 GHz) on 
the GBT. Observations to conduct searches for CO emission 
in high- z sub-mm galaxies, guided only by their photometric 
redshifts, are currently underway on the GBT. 



6 CONCLUSIONS 

We have derived meaningful photometric-redshift probabil- 
ity distributions for individual sub-mm galaxies using Monte 
Carlo simulations that take into account a range of measure- 
ment errors and model-dependent uncertainties. From the 
individual redshift probability distributions, we have mea- 
sured the cumulative redshift distribution for the sub-mm 
population of galaxies identified in blank field surveys. We 
find that the multi-wavelength FIR-radio data for 50 sub- 
mm galaxies detected (> 3.5<r) at 850/im are consistent with 
~ 50—90% of the population lying at redshifts z = 2—4, with 
only a small fraction, ~ 10% of the population at z < 2. We 
also show that up to 50% of the galaxies have colours that 
are consistent with z > 4. The possibility that the sub-mm 
galaxy population contains a significant fraction undergoing 
a major burst of star formation at such early epochs has im- 
portant consequences for models of galaxy formation. This 
suggestion is still highly dependent on the prior assumption 
of the evolutionary model for the sub-mm population which 
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affects most strongly those sub-mm sources detected only at 
850/wn. 

Shorter-wavelength sub-mm data (250 — 500/im) from a 
future balloon-borne experiment, BLAST, will provide pow- 
erful additional constraints (Az ~ ±0.5) on the redshift dis- 
tribution of all the SCUBA and MAMBO galaxies at z > 2. 
A photometric-redshift accuracy of ±0.5 for an individual 
galaxy will be sufficient to provide robust statistical mea- 
surements of the redshift distribution of the entire popula- 
tion, the global history of dust-obscured star formation and 
the clustering properties of sub-mm galaxies. 

We have argued that the large dispersion in the radio 
to sub-mm SEDs (or 850/im/1.4GHz colour) of the template 
galaxies naturally translates into an imprecise measure of 
the photometric redshift, and that radio observations can- 
not be used efficiently to measure the individual redshifts 
of thousands of galaxies detected in future wide-area sub- 
mm surveys. The sub-mm/radio spectral index is still able 
to discriminate between the z < 2 and z > 2 regime. 

We have an on-going programme to conduct centimetre- 
wavelength observations of CO J=l-0 molecular line tran- 
sitions in high- 2 sub-mm galaxies on the 100-m GBT. The 
observations, however, first require an accurate photometric- 
redshift probability distribution to guide the choice of re- 
ceiver and tuning parameters. We illustrate the power of 
our photometric-redshift technique with the example of 
HDF850.1, one of the best-studied sub-mm galaxies (Hughes 
ct al. 1998; Downes et al. 1999; Dunlop et al. 2002). 

Within the next few years, the combination of ground- 
based, balloon-borne and satellite experiments will provide 
accurately calibrated rest-frame sub-mm-FIR data for a sig- 
nificant population of sub-mm galaxies. The centimetre and 
millimetre wavelength spectroscopic redshift confirmation of 
their photometric redshifts will calibrate the method out- 
lined in this paper. Thus we can finally look forward to 
breaking the "redshift deadlock" that currently prevents an 
accurate understanding of the nature of the sub-mm galaxy 
population and their evolutionary history. 
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APPENDIX A: REDSHIFT ESTIMATES FOR 
INDIVIDUAL SUB-MM GALAXIES 

A catalogue of redshift probability distributions and SEDs 
for the galaxies detected in more than 2 bands (Tables 1 and 
2) is included. We also include the corresponding colour- 
colour-redshift or colour-flux-redshift plots from which they 
have been derived. All estimates are derived for model le2. 

For galaxies detected in just one band (Table 3) we 
present a selection of the results derived for 6 out of the 
46 galaxies included in Table 3. The general uniformity of 
their redshift distributions is a fair representation of those 
galaxies not included in this appendix. 
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